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ABSTRACT 

X-ray  diffraction  was  utilized  to  determine  the  root  mean  square  (r.m.s.)  strains  and 
average  particle  sizes  in  evaporated  gold  thin  films  on  aluminum  nitride  substrates  as  a 
function  of  substrate  surface  condition  prior  to  deposition.  The  substrate  treatments 
evaluated  were  surface  roughness,  use  of  titanium  and  chromium  inter-layers,  presence  of 
an  oxide  layer  on  the  substrate  surface  and  vacuum  conditions  used  during  deposition. 
The  Warren-Averbach  method  was  utilized  to  obtain  the  r.m.s.  strains  and  particle  sizes 
from  peak  breadth  data,  using  both  cosine  and  modulus  methods.  It  was  concluded  that 
the  highest  strain  deviations,  and  therefore,  the  largest  film  plastic  deformation,  occurred 
when  the  substrate  surface  was  rough,  when  chromium  was  used  as  an  inter-layer  and 
when  ultra  high  vacuum  conditions  were  used  during  deposition.  It  is  proposed  that  for  a 
fixed  film-substrate  system,  the  r.m.s.  strain,  which  is  indicative  of  the  level  of  plastic 
deformation  in  the  film  due  to  differential  contraction  following  deposition,  may  serve  as 
an  indirect  measure  of  the  interfacial  adhesion. 
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I.  INTRODUCTION 


Since  the  electronic  revolution  of  the  1960's  and  70's,  the  electronics  industry  has 
striven  for  ever  smaller  and  ever  more  powerful  electtonic  circuits.  Figure  1  charts  the 
evolution  of  circuit  power  dissipation  requirements  over  time.  The  steadily  increasing 
demand  on  small  circuit  size,  cost  and  reliability  results  in  the  constant  search  for 
improved  heat  dissipation  designs.  Since  an  integrated  circuit  must  typically  be  maintained 
below  80  to  100°C  for  component  reliability,  and  is  usually  subjected  to  numerous  thermal 
cycles  during  its  lifetime,  proper  design  of  packages  for  modern  electronic  components 
pose  some  serious  challenges. 

One  of  the  principal  properties  of  concern  is  the  thermal  conductivity  of  the  substrate 
on  which  chips  are  mounted.  A  higher  thermal  conductivity  allows  for  increased  heat 
removal  capability  and  a  higher  circuit  density.  Aluminum  nitride  (A1N)  offers  and 
attractive  choice  for  chip  substrates  over  other  substrates  such  as  alumina  (A1203)  and 
beryllium  oxide  (BeO).  Table  1  lists  engineering  properties  of  importance  of  materials 
currently  used  as  chip  substrates  and  how  they  compare  to  A1N.  Although  Table  1  lists 
the  thermal  conductivity  of  A1N  as  ranging  from  100-170  W/K-m,  it  has  a  theoretical 
value  of  over  300  W/K-m  and  other  attractive  properties  such  as  cost,  ease  of 
manufacture,  non  toxicity  and  good  insulating  capacity  [Ref.  3,  4],  which  make  it  very 
attractive  as  a  substrate  material  for  ceramic  packages. 

Metallization  (the  placement  of  a  thin  conductive  metal  film  on  the  substrate  surface) 
is  a  vital  step  in  providing  interconnects  between  devices,  and  for  the  attachment  of  seal 


rings  to  the  substrate.  The  metal  film,  which  may  be  applied  to  the  substrate  by  thin  film 
processes  (such  as  vacuum  evaporation  or  sputtering)  must  adhere  to  the  substrate  under 
the  presence  of  thermal  stresses  caused  by  thermal  cycling  without  cracking  or  peeling. 

Recent  attempts  to  improve  copper  and  gold  film  adhesion  to  A1N  substrates  have 
been  successful,  but  involve  application  of  oxide  layers  or  use  of  metal  pastes  that  typically 
lower  the  overall  thermal  conductivity  of  the  product.  Gold,  being  inert  and  possessing 
good  conductive  properties,  is  a  commonly  used  metallization,  however,  little  information 
is  available  on  the  influence  of  various  process  parameters  on  the  adhesion  of  gold  films  to 
A1N  substrates.  Especially  little  systematic  information  is  available  on  the  impact  of 
substrate  surface  conditions  prior  to  metallization  on  adhesion  and  the  state  of  stress  of 
the  film. 

Since  X-ray  diffraction  (XRD)  is  completely  nondestructive,  it  has  been  used 
extensively  to  study  strains  and/or  stresses  in  crystalline  thin  films  through  the  study  of 
diffraction  peak  shifts  and  peak  broadening.  While  techniques  based  on  peak  shifts  yield 
the  average  long-range  strains/stresses  in  materials  [Ref.  5],  those  based  on  peak  breadths 
yield  a  measure  of  strain/sttess  dispersion  (i.e.,  the  deviation  of  the  stains/stresses  from 
the  mean  long-range  value).  Of  the  techniques  based  on  line  breadth,  the  Fourier  method 
of  Warren  and  Averbach  [Ref.  6]  is  best  characterized  mathematically,  and  yields  a 
measure  of  both  the  average  crystallite  size  (which  may  represent  grain/subgrain  size  or 
inter-dislocation  spacing)  and  the  root  mean  square  (r.m.s.)  strain,  which  represents  local 
variations  of  the  lattice  strain  from  the  mean  long  range  strain  value.  Since  these 
deviations  usually  occur  due  to  plastic  deformation,  the  strain  dispersion  can  be  used  as  a 
measure  of  the  degree  of  plastic  deformation  of  the  films  due  to  differential  thermal 
contraction  of  the  film  and  the  substrate  during  cooling  following  deposition.  Further,  for 
a  given  pair  of  film  and  substrate  materials,  better  interfacially  bonded  systems  are 


expected  to  undergo  a  larger  degree  of  plastic  deformation,  therefore,  strain  dispersions 
obtained  from  peak  breadths  may  serve  as  an  indirect  measure  of  interfacial  adhesion. 

Accordingly,  the  goal  of  this  study  has  been  to  use  the  Warren  and  Averbach 
technique  to  determine  crystallite  size  and  r.m.s.  strains  following  metallization  using 
different  process  conditions  (primarily,  different  substrate  surface  conditions)  to 
qualitatively  examine  the  influence  of  these  conditions  on  metallization  stresses  and 
adhesion. 
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Figure  1.  Chip  Power  Dissipation  Over  Time  [Ref.  1] 


Table  1.  PROPERTIES  OF  ALN  AND  OTHER  CERAMIC  SUBSTRATES  [Ref.  2] 


AIN 

A1203 

BeO 

SiC 

Thermal  Cond. 
(W/m-K) 

100-170 

20 

250 

270 

Coeff.  Th.  Exp. 

(xio-Vo 

4.5 

7.3 

8.0 

3.7 

Dielectric  Str. 
(20°C)(kV/mm) 

14-17 

10 

10-14 

0.07 

Density  (g/cm3) 

3.3 

3.9 

2.9 

3.2 

Bend  Str. 
(MPa) 

3.9-4.9 

2.9 

1.6-2.3 

4.4 

Young  Mod. 
(xl07)(MPa) 

2.81 

3.7 

3.2 

4.0 

Sintering  Temp. 
(°C) 

1800 

1500 

2000 

II.  BACKGROUND 

A.     ALUMINUM  NITRIDE 

1.  Microstructure 

The  ceramic,  Aluminum  Nitride  (A1N)  is  a  covalent  compound  formed  by  the 
reaction  of  aluminum  and  nitrogen.  After  formation,  A1N  crystallizes  in  a  hexagonal 
wurtzite  structure  as  shown  in  Figure  2.  A1N  possesses  a  closed  packed  A1N  double  layer 
and  is  not  known  to  form  polytypes,  wherein  the  lattice  takes  on  a  cubic  structure  [Ref.  8]. 
The  microsttucture  can  be  severely  altered  in  the  presence  of  oxygen  whereby  the  oxygen 
replaces  three  nitrogen  atoms  at  the  exclusion  of  one  aluminum  atom  [Ref.  9].  The  effect 
of  oxygen  on  the  lattice  structure  is  of  importance  when  considering  metallization 
treatments  and  the  thermal  properties  of  A1N  and  is  addressed  later.  General  properties  of 
A1N  are  given  in  Table  1. 

2.  Fabrication 

There  are  several  techniques  available  to  form  A1N  powder  allowing  for  great 
flexibility.  Additionally,  the  raw  materials  needed  to  produce  A1N  powder  are  abundant 
and  relatively  inexpensive.  The  most  common  ways  to  commercially  make  A1N  powder  is 
by  the  reaction  of  ammonia  with  aluminum  chloride  or  by  direct  nitration  of  metallic 
aluminum  followed  by  carbon  reduction  in  a  nitrogen  environment  [Ref.  10].  More  recent 
methods  include  chemical  vapor  deposition,  polymer  pyrolysis  and  arc  plasma  techniques 
[Ref.  11,  12]. 

A  critical  requirement  in  the  production  of  useful  A1N  powder  is  the  purity  level 
achieved  during  production.  Excessive  amounts  of  carbon  or  oxygen  contaminants  can 
reduce  the  thermal  conductivity  of  A1N  product  following  sintering.  Additionally,  the 
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Ficure  2.  Lattice  Structure  for  A1N  [Ref.  7] 


presence  of  contaminants  can  hinder  the  adhesion  of  metal  thin  films  to  the  substrate 
surface.  The  latest  techniques  currently  in  use  have  resulted  in  A1N  powder  oxygen  levels 
as  low  as  1%  and  carbon  levels  as  low  as  0.02  weight  percent  [Ref.  13,  14]. 

The  method  chosen  to  manufacture  the  A1N  powder  also  affects  the  resulting 
particle  size,  a  key  parameter  in  the  sintering  process.  Current  techniques  offer  a  range  of 
grain  sizes  from  as  small  as  30  nm  up  to  8  microns  [Ref.  15]. 

Following  powder  production,  sintering  is  required  to  form  the  final  product. 
Several  sintering  aids  have  recently  been  studied  with  varying  effectiveness.  The  most 
common  sintering  aids  include  rare-earth  or  alkaline-earth  oxides,  such  as  calcium  oxide, 
yttrium  fluoride,  calcium  fluoride  and  yttria  [Ref.  16,  17].  Additionally,  carbon  doping  is 
often  performed  during  the  sintering  operation  [Ref.  18].  The  choice  of  sintering  aids  and 
doping  agents  have  a  significant  impact  on  oxygen  diffusion,  final  particle  size, 
densification  and  thermal  conductivity  [Ref.  19].  Some  recent  methods  have  resulted  in  a 
densification  as  high  as  99%  of  theoretical  and  a  thermal  conductivity  as  high  as  260  W/K- 
m,  a  nearly  ten  fold  increase  when  compared  to  the  thermal  conductivity  of  alumina  [Ref. 
20,21]. 

3.      Metallization 

Metallization  of  A1N  substrates  is  not  a  perfected  procedure  and  is  under  current 
study  by  a  variety  of  institutions  and  companies.  While  some  successes  have  been 
achieved  in  improving  the  adhesive  strength  of  metal  thin  films  to  the  substrate  surface, 
this  often  results  in  a  reduction  in  thermal  conductivity.  Processes  currently  under  study 
are  now  briefly  examined. 

Co-firing  A1N  with  tungsten  and  with  Mo-Mn  using  fritted  pastes  has  resulted  in 
good  adhesion  strengths,  relying  on  metal  diffusion  and  subsequent  mechanical 
interlocking  [Ref.  22].  Normal  adhesion  strengths  greater  than  20  MPa  have  been 


reported  using  this  method,  but  due  to  the  low  electrical  conductivity  of  W  and  Mo, 
resistive  heating  effects  have  reduced  the  thermal  conductivity  by  19%  and  9% 
respectively  [Ref.  23,24]. 

Thick  film  metallization  of  A1N  substrates,  using  commercially  available  pastes 
of  Au,  Cu  and  Ag-Pd,  has  similarly  proven  to  be  successful,  resulting  in  normal  adhesion 
strengths  from  about  14  MPa  to  18  MPa,  but  is  limited  by  the  formation  of  either  a  glassy 
phase  or  a  mixed  oxide  layer  near  the  interface  that  results  in  an  insulating  effect  [Ref.  25- 
29]. 

Due  to  the  absence  of  a  glassy  phase,  which  could  otherwise  be  exploited  for 
bonding,  metallization  using  fritted  pastes  of  various  metallic  oxides  and  glassy  pastes, 
have  been  performed  resulting  in  normal  adhesive  strengths  as  high  as  high  as  20MPa. 
However,  problems  encountered  with  these  methods  include  the  need  for  exceptionally 
pure  pastes,  chemical  and  thermomechanical  incompatibility  between  the  paste  and 
substrate,  oxygen  diffusion,  and  the  formation  of  a  glassy  inter-layer  which  causes  an 
insulating  effect  [Ref.  30].  Improper  control  of  the  paste  chemistry  can  result  in  poor 
adhesion,  and  presence  of  a  glassy  inter-layer  has  limited  the  achieved  thermal 
conductivity  to  about  150  W/K-m  using  these  methods  [Ref.  31]. 

Directly  bonded  copper  (DBC)  metallizations,  which  rely  on  Cu-0  bonding, 
have  been  successfully  applied  to  A1N  substrates  yielding  peel  strengths  as  high  as  8.5 
kg/cm  [Ref.  32].  The  success  of  the  DBC  method  was  found  to  be  strongly  dependent  on 
the  thickness  of  a  pre-applied  oxide  layer  that  promotes  Cu-0  bonding,  but  lowers  the 
thermal  conductivity  of  the  finished  product  due  to  the  formation  of  A1203  and  Cu20,  both 
of  which  have  a  thermal  conductivity  lower  than  that  of  A1N  [Ref.  33]. 

Oxygen  diffusion  can  be  a  problem  when  processing  A1N  for  metallization,  even 
at  near  ambient  temperatures,  since  the  surface  of  the  substrate  can  quickly  pick  up 


oxygen  from  the  air  forming  A1203  at  the  interface  in  a  matter  of  hours  [Ref.  34]. 
Impurity  pickup,  especially  carbon  pickup,  is  also  a  handling  concern.  Methods  to  reduce 
hydrocarbons  at  the  substrate  surface,  therefore,  is  an  important  processing  requirement. 

With  respect  to  the  methods  of  metallization  outlined  above,  metallization  of 
A1N  with  highly  conductive  films,  such  as  Au  or  Cu,  appear  to  hold  promise  in  both 
providing  good  electrical  conductivity  as  well  as  avoiding  a  reduction  in  the  thermal 
conductivity,  if  they  can  be  reliably  adhered  to  the  substtate  surface. 

Substrate  surface  roughness  is  likely  to  influence  mechanical  interlocking  at  the 
interface,  hence,  the  adhesion  of  thin  films  to  substrates.  However,  no  systematic  study  of 
the  effect  of  substrate  surface  finish  on  stress  and  adhesion  in  films  on  A1N  has  been 
reported  to  date.  Therefore,  substrate  surface  roughness  was  a  chosen  variable  in  this 
study.  Unoxidized  A1N  substrates  of  varying  surface  roughness  metallized  with  pure  gold 
were  examined  using  XRD. 

Studies  of  the  metal-AIN  interface  have  shown  that  the  interfacial  bond  is 
controlled  by  a  metal-AIN  interaction  vice  a  metal-0  interaction  in  A1203  metallizations 
[Ref.  35].  The  metal-AIN  interaction,  and  thus  the  interfacial  adhesion,  can  often  be 
strengthened  by  depositing  a  thin  inter-layer  of  a  reactive  metal  (e.g.,  Ti)  between  the 
metallization  and  the  A1N  substrate  [Ref.  36-38].  Ti  has  shown  to  enhance  wetting  at 
metal-AIN  interfaces  via  formation  of  TiN  and  metallic  aluminum.  Because  of  these 
advantages,  the  TiN  system  has  been  the  focus  of  most  studies  of  thin  film  metallization  to 
date,  yet  little  quantitative  data  on  the  effect  of  Ti  on  the  strength  of  metal-AIN  interfaces 
is  available.  Large  residual  stresses  have  been  known  to  be  present  in  thin  metal  films  on 
ceramic  substtates  for  a  long  time  (with  resultant  implications  on  the  mechanical 
properties  and  the  reliability  of  the  interface)  [Ref.  39],  however,  no  systematic  study  of 


the  dependence  of  interfacial  adhesion  on  the  stresses  at  thin  film/AIN  interfaces  has  been 
reported. 

In  addition  to  a  Ti  inter-layer,  a  Cr  inter-layer  will  be  investigated.  Unlike  Ti,  Cr 
does  not  form  any  intermetallics  with  Au  yet,  Cr  forms  a  solid  solution  with  Al  and  several 
CrmAln-type  intermetallics  when  A1N  and  Cr  are  heated  to  1000°C  [Ref.  40,  41].  This  can 
potentially  give  rise  to  good  bonding  between  Cr  and  A1N,  due  to  the  absence  of 
intermetallics  in  the  Au-Cr  system,  and  may  yield  a  metallization-substrate  interface  with 
better  interfacial  mechanical  properties  than  the  AIN-Ti  system.  Therefore,  A1N  substrates 
metallized  with  Au,  using  Cr  and  Ti  as  inter-layers,  are  investigated  in  this  study. 

As  outlined  above,  formation  of  an  oxide  layer  has  been  shown  to  improve 

bonding  in  Cu  thin  films  applied  to  A1N  substrates  through  Cu-O  bonding.  It  is  known 
that  Au  forms  a  metal  oxide,  Au~02,  at  temperatures  above  500°C  [Ref.  42].  However, 

the  effect  of  the  presence  of  a  thin  layer  of  A1203  on  the  A1N  surface  on  bonding  with  Au 
metallization  is  not  known,  and  is  investigated  here. 

Thin  film  metallizations  are  usually  deposited  on  substrates  with  ultra  high 
vacuum  (UHV)  or  high  vacuum  (HV)  conditions.  Generally,  higher  vacuum  promotes 
better  epitaxial  (oriented)  crystal  growth  with  the  vacuum  playing  a  key  role  in  the  thin 
film  microstructure  [Ref.  43].  Presently,  little  information  on  the  effect  of  metallization 
pressure  on  the  film  characteristics  when  direcdy  applying  Au  to  A1N  substrates  is 
available.  Therefore,  metallizations  under  UHV  and  HV  conditions  were  performed. 

Finally,  substrate  temperature  during  metallization  is  also  known  to  profoundly 
affect  epitaxial  growth  of  thin  films  and  may  also  play  a  key  role  in  reaction  metal  bonding 
at  the  interface  [Ref.  44].  Therefore,  metallizations  at  two  different  substrate 
temperatures  using  unoxidized  A1N  substrates  with  Cr  inter-layers  was  performed. 
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B.     STRAIN/CRYSTALLITE  SIZE  MEASUREMENT  BY  X-RAY 
DIFFRACTION:  SUMMARY  OF  THE  WARREN-AVERBACH 
METHOD 

1.      Representation  of  a  Periodic  Function  by  a  Fourier  Series 

A  complete  description  and  derivation  of  the  Warren-Averbach  method  is 
available  in  [Ref.  6,  45,  47].  The  summary  of  the  Warren  Averbach  method  presented  in 
the  following  sections  is  based  on  [Ref.  45,  46]. 

The  first  step  in  obtaining  quantitative  results  from  X-ray  diffraction  by  the 
Warren-Averbach  method  requires  representing  the  diffraction  profile  by  a  Fourier  series. 
Since  any  periodic  function  can  be  represented  by  a  Fourier  series,  one  can  assume  a 
measured  diffraction  profile,  designated  /(*),  to  be  a  section  of  a  periodic  function  having 
the  shape  of  the  measured  profile  as  shown  in  Figure  3.  It  can  be  shown  that  this  function 
can  be  represented  by 


/(*)=  ^    \cos2nn  —  \  + Bns\n2nn  — 


(1) 


where  a  is  the  period  of  the  function,  —  to  —  is  the  interval  of  interest  and  n  takes  on 

integer  values  and  is  referred  to  as  the  harmonic  number.  It  can  also  be  shown  that  the 
Fourier  coefficients  An  and  Dn,  can  be  solved  for  by 


+— 

2 


An=-\f(x)cos2Kn(-\h  (2) 


+— 
Bn=jjf(x)sm2nn(±\k.  (3) 
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Figure  3.  Representation  of  a  Periodic  Function  [Ref.  47] 


12 


Alternatively,  equation  (1)  can  be  expressed  in  its  complex  form  as 

/(*)=  fjC(n)e-2mKxla)  (4) 

n=-<»> 

where  the  complex  coefficients,  C„,  are  composed  of  a  real  and  an  imaginary  part  such 
that 

C{n)  =  Cr(n)  +  iC,(n) 
where  the  subscript  r  refers  to  the  real  part,  and  the  subscript  i  refers  to  the  imaginary 
part. 

Since  the  term  e~2w"(Jt/a)  in  equation  (4)  can  be  expressed  in  its  trigonometric 
form  as 


-2rcm  -  (  x\  (  X 

e       Ka' =cos2nn\  —  \-isin2nn\  —  \, 
then 


/(*)=  ]T    Cf(AOcos27i«f-)  +  CK/Osin27r/2f-l  +  /C,(/Ocos2roif-l-/Cr(«)sin2K«f-l 


•  (5) 


In  equation  (5),  the  imaginary  terms  must  cancel  since  f(x)  is  real.  This  can 
only  happen  if  C(-n)  =  C*(n),  where  C*(n)  represents  the  complex  conjugate  of  C(n). 
This  result  allows  Cr(n)  and  Q(n)  to  be  expressed  as 


C,(«)  =  -  \f(x)cos2nn\-)dx  (6) 

a  J  V  a  J 
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and 


Ci(n)--\f(x)sm2nn(-\bc.  (7) 


Comparing  equations  (6)  and  (7)  with  equations  (2)  and  (3),  it  is  evident  that 
Cr(n)  =  An       and      C,(n)  =  Bn 
where  Cr(n)  represents  the  real  Fourier  coefficients  and  C,(/i)  represents  the  complex 
Fourier  coefficients.  Thus,  a  measured  diffraction  profile  is  can  be  expressed  by  a  Fourier 
series  with  real  and  complex  coefficients. 

2.      The  Measured,  Instrumental  and  Diffraction  Profiles 

The  measured  peak  profile,  designated  h(x),  comprises  the  insu'umental  profile 
in  conjunction  with  the  true  diffraction  profile.  The  instrumental  broadening  effect  is 
caused  by  incident  radiation  that  is  not  truly  monochromatic  and  the  fact  that  the  incident 
radiation  is  spread  across  a  finite  width  of  the  specimen,  vice  being  confined  to  an 
infinitely  thin  line  [Ref.  48].  This  instrumental  effect  is  apparent  even  when  the  specimen 
is  strain  free  and  of  a  sufficient  particle  size  that  avoids  particle  size  broadening.  The 
measured  profile  then,  is  a  convolution  of  the  instrumental  profile,  designated  g(x),  and 
the  pure  profile,  designated  f(x).  Figure  4  represents  the  situation  graphically.  Obtaining 
the  instrumental  profile,  g(x),  requires  obtaining  a  diffraction  profile  of  a  strain  free 
standard  of  a  sufficiently  large  particle  size,  so  that  neither  particle  size  broadening,  nor 
strain  broadening  is  apparent.  The  combined  effect  of  g(x)  and  f(x)  can  be  represented 
by  a  convolution  integral 

h(x)  =  -\g{z)f{x-z)dz  (8) 

A  J 
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Figure  4. (a)  Plot  of  Measured  Profile,  h(x) 


500 


Figure  4.(b)  Plot  of  Instrumental  Profile,  g(z) 


800  T 


Figure  4.(c)  Plot  of  True  Profile,  f(y) 
Figure  4.  Graphical  Representation  of  a  Deconvoluted  Profile  [Ref.  49] 
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where  A  is  the  area  under  the  curve  f(y). 

If  convenient  shapes  (or  functions)  of  the  three  profiles  are  assumed  (e.g. 
Cauchy  or  Gaussian),  the  relation  between  the  integral  line  breadths,  B,  of  the  three 
profiles  can  be  expressed  as 

Cauchy:  u  =  U/(\  +  k2u2) 


and 


B2(h)  =  B2(g)  +  B2(f) 


Gaussian:  u  =  Uek~" 


B(h)  =  B(g)+B(f)  [Ref.  50]. 
Actual  diffraction  profiles,  however,  are  neither  purely  Gaussian  nor  purely 
Cauchy,  making  a  Fourier  representation  of  the  profiles  and  subsequent  deconvolution 
necessary. 

3.      The  Fourier  Transform  or  Stokes  Deconvolution 
Representing  h(x),  f(y)  and  g(z)  in  Fourier  series  form 


-27Cn; 


f(y)  =  JJF(n)e      {a)  (9.1) 


^)  =  XG(n>       {a)  (9.2) 


h(x)  =  ^H(n")e        {a)  (9.3) 


where:  F,  G  and  H  are  the  Fourier  coefficients  and  n ,  n'  and  n"  are  the  harmonic 
numbers.  Substituting  equations  9.1,  9.2  and  9.3  into  equation  8  yields 
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h(x)  =  —  J  £G(n>       u;^F(«)e       v  a  >dz 


where  the  limits  of  integration  can  be  replaced  by  —  to  —  so  that  g(z)  is  non-vanishing 

only  in  this  interval.  Combining  all  terms  containing  z,  and  replacing  the  limits  of 
integration  yields 


1  -Jn/i)    f    -2*i(«'-/0f-l  , 

/j(x)  =  — XXG(n')F(«)e       UjJe  u'<fe 

n'      n  a 

Since, 


-2jB(n'-/i)f-]         f  a  when  n'  =  n 
0  when  n'  *  « 


r2    -2jb(»'-«)T— ]  f 

a  t 

~2 

and  the  condition  of  interest  is  ri  =  n,  then 


h(x)  =  —  YJG(n)F(n)e       Kai      (for/i  =  0  to~).  (10) 


Replacing  /i(jc)  in  equation  (10)  by  equation  (9.3),  yields 


2^//(n)e       u;  =— 2^G(/i)F(«)e       u; 


or 


//(«)  =  4g(«)F(«).  (11) 

A 
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Since  only  profile  shape,  not  profile  size,  is  of  interest,  the  constant  —  can  be  dropped 

A 


leaving 


H{n)  =  G{n)F{n). 


(12) 


This  relation  can  be  rewritten  in  complex  form  as 


F,(/t) +  //?(«)  = 


Hr(n)  +  iH(n)i 


Gr(n)  +  iGt(n) 

Multiplying  the  numerator  and  denominator  by  the  complex  conjugate  of  G(n)  and 
equating  the  real  and  imaginary  parts  yields 


rAn)  = : z 

G>)  +  G,2(n) 


(13.1) 


and 


Ft(n)  = 


//,-(/OGr(«)-//r(;QG,.(fO 


(13.2) 


From  the  two  experimentally  measured  curves  h(x)  and  g(z),  the  real  (cosine) 
coefficients  (Hr(n)  and  Gr(n))  and  the  imaginary  (sine)  coefficients  (//,-(«)  and  G-(«))  can 
be  determined.  The  broadened  diffraction  profile  cosine  and  sine  coefficients  are  then 
determined  using  equations  (13.1)  and  (13.2)  allowing  the  pure  diffraction  profile  to  be 
constructed  from 


/()0  =  X  ^V('0cos2tw  —  +  F/(«)sin27t«  — 


(14) 
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4.      Determination  of  Lattice  Strains  and  Crystallite  Size 

Cold  worked  metals  may  contain  numerous,  non-random  dislocation  arrays, 
such  as  low  angle  grain  boundaries  which  subdivide  the  original  grains  into  smaller 
coherent  domains  each  slightly  misoriented  from  the  next.  Therefore,  non-vanishing 
diffraction  intensities  can  occur  at  angles  other  than  that  satisfying  the  exact  Bragg 
condition,  resulting  in  broadening  of  the  diffraction  line.  Such  peak  broadening  can  also 
occur  if  the  crystal  is  divided  into  smaller,  relatively  misoriented  crystallites  by  random 
dislocation  arrays,  or  if  the  crystal  particle  sizes  are  very  small.  Typically,  crystallite  sizes 
of  100  - 1000  A  yield  significant  diffraction  broadening  [Ref.  51].  If  the  broadening  is  due 
only  to  this  crystallite  size  effect,  the  average  crystallite  size  is  given  by 


D-F^  (15) 

pcos9 


where:   Dv  is  the  crystallite  size  in  angstroms,  k  is  a  proportionality  constant  close  to  1.0, 

P  is  the  integral  breadth  of  the  true  profile  such  that  p  = — ,  9  is  the  Bragg 

peak  maximum 

angle  and  X  is  the  wavelength  of  the  incident  x-rays  in  A  [Ref.  52].  The  dislocations, 

which  produce  the  subdivisions  into  domains,  can  also  produce  both  tensile  and 

compressive  strains  within  each  domain.  These  tensile  and  compressive  strains  result  in  a 

deviation  of  the  local  lattice  spacings  from  the  mean  lattice  spacing,  thus,  resulting  in  peak 

broadening. 

If  the  broadening  is  due  to  microstrains  only,  the  mean  lattice  strain  is  given  by 


e  =  — L  [Ref.  53].  (16) 

4  tan  9 
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Usually,  however,  both  crystallite  size  and  strain  broadening  effects  are  convoluted  in  the 
diffraction  profile,  and  a  means  of  deconvolution  is  required  to  determine  the  respe  tive 
contributions  to  the  integral  breadth.  Because  of  the  difference  in  functionalities  of  Dv  and 
e,  with  respect  to  9,  their  contributions  to  (3  can  be  deconvoluted.  From  equation  (14), 


/O0  =  X  Fr(n)cos2wif^+/?(ii)sin2jwKJ 


(17) 


Utilizing  a  Stokes  deconvolution  of  the  measured  profile,  the  real  and  complex  Fourier 
coefficients,  of  the  diffraction  profile  Fr(n)  and  F^n),  can  be  obtained  using  equation  (13). 
Renaming  Fr(n)  as  An  and F^n)  as  Dn,  Warren  and  Averbach  were  able  to  determine  the 
complex  coefficients  as 


4=^-(cos27i/Z„>  (18) 


and 


Bn  =  -^(sm2nlZn)  (19) 


where: 


_  N,,  (sample)  _  no.  of  cells  in  the  whole  sample  having  an  na'  neighbor  in  the  same  column 
iV(column)  no.  of  columns  in  the  sample 

(which  can  be  rewritten  as  Nn  =  the  number  of  cells  with  an  nUl  neighbor  in  the  same 

column),  jV3  =  the  number  of  cells  per  column,  /*  =  the  harmonic  number,  /  =  the  order  of 

reflection,  Oj  =  the  unit  lattice  spacing  perpendicular  to  the  00/  reflections,  L  =  nai  =  the 
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real  space  distance  between  a  pair  of  n^  neighbors  along  a  line  perpendicular  to  the  00/ 
planes,  Z„  =  the  number  (or  fraction)  of  lattice  spacings  perpendicular  to  the  00/  reflection 
by  which  the  distance  between  a  pair  of  n^  neighbors  has  changed  due  to  distortion, 

AL  =  Znflj  =  the  change  in  distance  between  a  pair  of  nUl  neighbors  due  to  distortion.  With 

this  result,  the  strain,  eL,  can  be  expressed  as  zL  =  —  =  — --1 

L       n 

The  parenthesis  in  equation  (18)  and  equation  (19)  signify  average  values  since  no 

assumptions  about  the  direction  of  relative  displacements  of  the  unit  cells  is  made  in  the 

derivation.  The  cosine  coefficient  of  equation  (18)  is  the  product  of  a  size  related  term 

and  a  distortion  related  term  as 


AT=^-  (20.1) 

/V3 


and 

A?aonion=(cos2nlZn)  (20.2) 


the  former  being  dependent  on  the  column  length  only  (representing  the  size  component), 
and  the  latter  being  dependent  on  the  strain  only  (representing  the  distortion  component). 
The  cosine  coefficient  measured  can  now  be  represented  by  the  product 

An=A^A^stor,ion.  (21) 


'The  derivation  of  the  Warren  and  Averbach  analysis  is  begun  from  a  general  assumption  that  the  sample 
is  orthorhombic,  with  reflection  from  the  (00/)  plane.  It  is  assumed  the  sample  is  comprised  of  unit  cells 
arranged  in  columns,  with  spacing  L  =  nai,  between  an  arbitrary  pair  of  unit  cells  located  in  a  column 
whose  axis  is  normal  to  the  reflecting  plane  (see  Figures  5  and  6).  However,  the  analysis  is  completely 
general,  and  is  applicable  to  all  crystal  systems  and  reflections  from  non-(00/)  type  planes  [Ref.  54]. 
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Figure  5.  Real  Space  Representation  of  a  Columnar  Material 
used  in  the  Warren  Averbach  Derivation  [Ref.  55] 


Figure  6.  Reciprocal  Space  Representation  of  a  00/  Reflecting  Plane 
used  in  the  Warren  Averbach  Derivation  [Ref.  56] 
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Assuming  multiple  order  reflections  have  been  measured,  (e.g.  001,  002,  003), 
and  their  respective  Fourier  coefficients  have  been  determined,  it  is  possible  to  separate 
the  size  coefficients,  A^'ze ,  from  the  distortion  coefficients,  Afs,or"on ,  since  the  size 
coefficient  is  independent  of  /  ,  whereas  the  distortion  coefficients  depends  on  /  as  well 
as  the  relative  unit  cell  displacement  Z„.  Assuming  the  product  lZn  is  small, 
(cos27r/Z„)  — >  1  -  2k212(zI^  which  can  be  expressed  in  logarithmic  form  as 

ln(cos27i/Zn)  =  ln(l-27i2/2(z„2))  =  -27r2/2(Z;).  (22) 

Expressing  equation  (20)  in  its  logarithmic  form  with  substitution  from  equation 

(22),  yields 

lnA,(/)  =  ln^-27:2/2(z2).  (23) 

Plotting  \nAn(l)  versus  /2,  gives  a  linear  relation  where  the  size  coefficient,  A^ze  is 
determined  from  the  y-intercept,  and  the  slope  yields  the  value  for  -2jt2(z2).  Figure  7 

shows  an  example  of  such  a  plot.  Since  the  change  in  displacement  in  the  direction  normal 
to  the  reflecting  planes  is  given  by  AL  =  ^(z,,),  and  L  =  na^  (refer  to  Figure  3),  then 

the  strain  can  be  expressed  as  (zL)  =  —  =  -L-£^,  and  (Z„)  =  U2L)n2 .  Since  n  and  /   are 
known,  the  root  mean  squared  (r.m.s.)  sttain  (e2- )      may  be  determined  from  the  slope  of 


the  plots  over  several  (at  least  two)  orders  of  reflection  and  represents  the  r.m.s.  strain  in 
the  direction  normal  to  the  reflecting  plane.  It  must  be  remembered  that  the  r.m.s.  value 
so  measured,  actually  represents  the  difference  between  the  total  residual  strain  at  the 
microlevel  and  the  long  range  residual  strain.  Thus,  the  calculated  r.m.s.  strain  is  a 
measure  of  the  deviations  of  the  local  microstrain  from  the  average  long  range  residual 
strain  in  the  sample.  Mathematically  this  is  given  by 
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\Ei./meas  ~Ueiot        £longrange/     )  -  \(£tot  )    /       2\£toteiong  range/  +  \(£long  range)    /•  (2<*) 

Now, 

/p  res  V  res  /        \ 

\c-tot  /        elong  range    'X'-'L/' 

but,  (el)  =  0,  since  if  the  positive  and  negative  components  of  eL  do  not  cancel  one 
another,  they  will  merely  contribute  to  e[^ngrange  and  result  in  a  peak  shift,  but  no  peak 

broadening.  Therefore, 

-2fres        /Fres\  =  -2f£res        ^2 

long  range  \    tot  /  ^     long  range' 

and 

(e0nieaS  =((C)2)-(C„gra„ge)2  [Ref-  58].  (25) 

Thus,  it  is  important  to  remember  that  the  r.m.s.  value  of  eL, which  is  obtained  by  the 
Warren- A verbach  method,  determines  not  the  average  strain,  but  rather,  the  deviation 
from  the  mean  (long  range)  sttain. 

The  method  described  achieves  the  desired  separation  using  the  real  (cosine) 
Fourier  coefficient,  An,  only,  and  is  the  traditional  cosine  method  of  arriving  at  the  r.m.s. 
strain  value.  This  method  assumes  that  the  origin  of  the  Fourier  series  expansion  is  the 
centroid  of  the  standard  and  measured  profiles,  and  that  each  centroid  location  is  known 
with  high  accuracy.  If  this  is  the  case,  the  imaginary  parts  of  the  Fourier  coefficients  are 
small,  resulting  in  reliable  results.  In  practice,  small  deviations  from  the  true  centroid 
positions  of  the  measured  and  standard  profiles  are  likely  to  result  (due  to  interference 
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between  peaks  or  due  to  background  errors)  causing  significant  error  to  the  calculated 
Fourier  coefficients  for  the  [rut  (deconvoluted)  profile  [Ref.  59]. 

A  second  method,  known  as  the  modulus  method,  uses  the  modulus  of  the  real 
and  imaginary  Fourier  coefficients  of  the  measured  and  structural  profiles  which  make  the 
location  of  the  origin  of  the  Fourier  series  expansion  irrelevant.  The  discussion  of  the 
modulus  method  that  follows  is  borrowed  from  [Ref.  60].  The  modulus  from  equation 
(13)  is  defined  as 


1/2 

\M(n)\= 


(Fr(nf  +  FM)2 


The  derivation  from  this  point  is  carried  out  in  the  same  manner  as  the  cosine  method 
resulting  in  the  relation 

\M(n)\=Kze -K«2n2l2{(zl)-(Zn)2). 

A  plot  of  \M(n)\  versus  I2  yields  A^,ze  as  in  the  cosine  method,  and  the  value 
{(z„2)  ~{Zn)  }m  which  can  be  expressed,  as  in  the  case  of  the  cosine  method,  as 
"l(ei)_(e)Ij     •  Since  (eL)  =  0  as  before, 

{((«=),)-(«=.^),r-(e1J)I 

which  represents  the  r.m.s.  strain.  The  particle  size  is  determined  in  the  same  manner  as 
the  cosine  method  discussed  below. 
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For  the  determination  of  crystallite  size,  the  values  of  Asn,ze  as  a  function  of 
harmonic  number,  n,  are  read  off  from  the  plot  of  lnAn  versus  I2  (Figure  7).  Then  A"v  is 
plotted  against  n ,  as  shown  in  Figure  8.  It  can  be  shown  that  when  n  — >  0, 


(  dAtv  > 


dn 


Ai->0  ^3 


(26) 


Therefore,  the  slope  of  the  Asn'ze  versus  n  plot  at  n  =  0  equals ,  N3  being  the  n  axis 

intercept  of  the  slope  of  the  Asn'ze  versus  n  curve  at  n  — >  0.  The  average  column  length 
perpendicular  to  the  00/  planes  (i.e.,  the  crystallite  size)  is  then  given  by  A^Oj,  where  <%  is 
the  unit  displacement  along  the  column  length. 

The  Warren- Averbach  approach  can  be  generalized  to  any  order  of  reflection  by 
powder  diffraction.  For  cubic  systems  equation  (23)  can  be  expressed  as 

2n2U2)L2(h2+k2+l2) 

\nAL(h0)  =  \nAr ^^— -2 (27) 

a 

where:  d  is  the  interplanar  spacing,  L  is  a  real  space  distance  along  the  columns  of  cells 

perpendicular  to  the  reflecting  planes,  and  a  is  the  lattice  constant.  A  graph  of  this 

relation  is  shown  in  Figure  9,  where  the  abscissa  intercept  of  the  initial  slope  represents  the 

average  column  length,  (D),  perpendicular  to  the  reflecting  planes. 

Observation  of  actual  A^  versus  n  plots  sometimes  show  a  concave  downward 

line  near  n  =  0,  which  has  come  to  be  known  as  the  "hook  effect"  [Ref.  63].  Since  the 

d2As 
second  derivative  of  Asn  with  respect  to  n  (i.e.,  — f-)  can  only  be  positive,  the  plot  of  Asn 

dn 
versus  n  can  be  concave  upward  but  never  concave  downward.  This  situation  is 
graphically  represented  in  Figure  9  by  the  dashed  line.  The  source  of  this  problem  can  be 
traced  to  using  line  profile  backgrounds  in  the  Fourier  analysis  that  are  too  high,  leading  to 
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Figure  7.  Plot  of  In  A^(l)    vs.  I2  [Ref.  57] 


Figure  8.  Plot  of  A£  vs.  n  [Ref.  61] 


<D> 

Figure  9.  Plot  of  A[  vs.  L  [Ref.  62] 
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unintended  peak  interference.  There  are  many  graphical  schemes  to  work  around  this 
problem,  but  so  far,  no  formal  way  has  been  developed.  The  easiest  way,  is  to  merely 
extrapolate  the  near  linear  portion  of  the  curve  back  to  the  ordinate  and  base  calculations 
on  this  value.  If  a  graphical  option  is  not  used,  then  background  adjustment  is  the  only 
other  way  to  avoid  this  effect. 

The  Warren-Averbach  method,  therefore,  allows  the  determination  of  particle 
size  and  strain  contribution  to  line  broadening  independently  and  simultaneously.  Though 
the  manual  calculations  are  laborious,  with  the  use  of  computer  driven  software  currently 
available  on  most  diffractometer  systems,  these  calculations  can  be  done  in  a  matter  of 
minutes  requiring  only  the  input  of  the  standard  diffraction  profile  and  the  measured 
diffraction  profile. 

In  the  above  analysis,  it  was  assumed  that  peak  broadening  occurs  due  to 
particle  size  and  strain  effects  only,  ignoring  any  contributions  from  deformation  or  twin 
faults. 
C.     PHYSICAL  SIGNIFICANCE  OF  THE  RESULTS  OBTAINED  FROM 

THE  WARREN-AVERBACH  ANALYSIS 

As  discussed  previously,  small  particle  size  can  contribute  to  line  broadening.  If  the 
sample  was  completely  free  from  dislocations  and  stacking  faults,  the  true  particle  size 
could  be  determined  directly  from  a  broadened  profile.  In  reality,  dislocations  can  be 
present,  to  a  high  degree,  in  random  or  non  random  configurations.  If  the  dislocations  are 
known  to  be  dispersed  randomly  throughout  the  material,  a  measure  of  the  particle  size 
would  be  an  indication  of  the  dislocation  density  of  the  material.  However,  dislocations 
that  are  not  randomly  dispersed  in  grains  and  form  low  angle  grain  boundaries  further 
subdivide  the  grain  into  smaller  subgrain  regions  leading  to  a  lower  realized  particle  size 
and  further  line  broadening.  Additionally,  dislocation  lines  break  up  the  ordered  planar 
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structure  causing  slight  disorientation  between  planes  resulting  in  a  range  of  angles 
satisfying  the  Bragg  condition  which  causes  line  broadening.  If  dislocations  are  known  to 
be  distributed  in  non  random  fashion  so  described,  the  particle  size  obtained  would  be 
representative  of  the  subgrain  size.  The  particle  size  determined  from  the  Warren- 
Averbach  method  must  therefore  be  corroborated  with  transmission  electron  microscopy 
(TEM)  studies  of  the  actual  film  microstructure.  TEM  corroboration  of  the  density  of 
stacking  faults  and  twins  should  likewise  be  conducted. 

Assuming  that  a  film  is  free  of  dislocations,  an  applied  stress  would  cause  a  uniform 
change  in  planar  spacing,  resulting  in  a  peak  shift.  The  total  strain  in  the  sample  would 
then  be  the  long  range  strain,  from  which  the  residual  stress  in  the  film  could  be  calculated. 
In  the  presence  of  dislocations  dispersed  in  the  material,  however,  local  strain  fields  may 
deviate  from  the  long  range  value.  Usually,  both  positive  and  negative  stress  states  will  be 
found  locally  due  to  the  presence  of  dislocations.  Even  if  the  average  of  these  stresses 
equals  zero,  and  no  peak  shift  is  noted,  the  r.m.s.  value  of  strain  as  measured  by  the 
Warren-Averbach  method  may  be  non-zero.  The  Warren  Averbach  analysis  yields  a  value 
of  strain  which  represents  the  difference  between  the  sum  of  the  positive  and  negative 
microstrains  squared  minus  the  long  range  (or  average)  sttain  squared.  This  result  then 
shows  the  possible  microsttain  deviation  from  the  average  strain,  not  the  average  strain 
itself.  The  stress  calculated  from  this  result,  therefore  represents  the  possible  sttess 
deviation  from  the  residual  stress  present  in  the  film.  Although  it  is  not  a  measure  of  the 
true  residual  stress,  it  is  none  the  less  useful  as  a  measure  of  the  degree  of  plastic 
deformation  of  the  film,  which  may,  under  certain  circumstances,  yield  a  measure  of 
interfacial  adhesion. 

Consider  two  film-substrate  pairs,  both  Au-AIN,  the  first  with  better  interfacial 
adhesion  than  the  second  due  to  differences  in  substrate  surface  treatment  prior  to 
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metallization.  Although  the  product  (ACTE)(AT)  for  the  given  materials  system  is  the 
same  for  both  pairs,  differences  in  the  strains  induced  in  the  films  are  expected  becc  jse  of 
different  degrees  of  strain  accommodation  at  the  interfaces.  If  plastic  deformation  occurs 
in  the  films  during  cooling,  the  stresses  induced  can  be  relaxed  to  a  level  below  the  yield 
strength  by  dislocation  generation.  The  residual  strain/stress  measured  via  peak  shift 
based  techniques  represents  the  relaxed  value.  If  both  samples  underwent  the  same  extent 
of  plastic  deformation,  then  the  su'ess  measured  by  peak  shifts  would  be  indicative  of  the 
interf acial  adhesion.  However,  in  general,  the  extent  of  plastic  deformation  can  be 
substantially  different  between  the  two  samples,  and  no  direct  correlation  between  the 
interfacial  adhesion  and  the  long  range  residual  stress/strain  measured  by  peak  shift  is 
obtainable. 

The  r.m.s.  strain  measured  based  on  peak  breadth,  on  the  other  hand,  represents  the 
deviation  of  the  local  microstrains  from  the  long  range  value,  and  thus,  is  indicative  of  the 
extent  of  plastic  deformation  of  the  film.  For  a  constant  (ACTE)(AT),  the  extent  of  plastic 
deformation  is  a  function  only  of  the  interfacial  adhesion  between  the  film-substrate  pair. 
Therefore,  the  r.m.s.  strain  is  likely  to  provide  an  indirect  measure  of  the  quality  of 
interfacial  adhesion.  The  condition  for  development  of  a  thin  film  metallization  procedure 
would  be  to  produce  a  film  with  good  adhesion  and  a  low  value  of  the  mean  long  range 
strain.  Having  a  low  strain  would  prevent  cracking  of  the  film  during  use,  whereas,  good 
interfacial  adhesion  would  prevent  peeling  of  the  film  at  its  interface.  The  existence  of  a 
low  strain  condition,  however,  does  not  necessarily  go  hand-in-hand  with  good  adhesion, 
requiring  both  the  long  range  strain  and  adhesion  to  be  evaluated  independently.  The  long 
range  stress/strain  can  be  evaluated  non-destructively  using  techniques  such  as  "the 
sin2  \j/"  measurement  method  [Ref.  64].  However,  nondestructive  evaluation  of  interfacial 

adhesion  is  not  feasible  at  present.  It  is  proposed  that  while  adhesive  properties  cannot  be 
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measured  directly  with  any  presently  available  x-ray  based  technique,  the  Warren- 
Averbach  method  may  provide  an  indirect  evaluation  of  the  extent  of  interfacial  adhesion, 
and  therefore  is  a  useful  tool  in  studying  the  quality  of  metallizations. 
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III.  EXPERIMENTAL 

A.  MATERIALS 

The  A1N  substrates  were  produced  and  provided  by  W.R.  Grace  and  Company  at  a 
purity  level  of  99.9%.  The  A1N  was  sintered  using  CaO  as  a  sintering  aid.  An  XRD  scan 
of  the  substrates  revealed  no  Ca-aluminate  peaks  that  would  otherwise  interfere  with  the 
Au  (1 1 1)  diffraction  line  after  metallization,  suggesting  that  the  CaO  resulting  from  the 
sintering  aid  additions,  was  lost  during  processing. 

The  Au  used  for  metallization  was  obtained  from  The  Aldrich  Chemical  Company 
and  was  used  in  the  powder  form  for  metallizations.  The  Au  powder  had  a  purity  level  of 
99.99%  and  an  average  particle  size  of  <44  microns  (-325  mesh).  This  powder  was  also 
used  to  prepare  the  Au  standard,  since  the  powder  size  was  sufficiently  large  (greater  than 
1  micron)  so  as  to  prevent  particle  size  broadening  of  the  standard  peaks  [Ref.  65]. 

B.  EQUIPMENT 

1.      Evaporator  Chamber  and  Sub  Components 

A  schematic  of  the  thin  film  evaporation  system  is  shown  in  Figure  10.  The 
chamber  was  made  for  metallizations  under  HV  conditions.  Six  vacuum  service  ports  with 
o-ring  seals  allow  for  a  variety  of  peripheral  equipment  and  monitoring  devices  to  be  used 
during  metallization.  Of  the  three  bottom  ports,  one  was  used  merely  as  a  viewing 
window  and  the  two  remaining  ports  were  used  for  power  feed-throughs.  An  external 
power  supply  was  used  to  supply  up  to  40  AC  amps,  single  phase,  current  to  the  power 
feed-throughs.  The  power  supply  was  equipped  with  an  output  amp  meter  and  rheostat 
that  allowed  for  the  precise  control  of  current.  This  allowed  the  operator  to  control  the 
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Figure  10.  Chamber  Layout  and  Equipment  Block  Diagram 
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metal  deposition  rate.  Heavy  gage  leads  with  friction  type  disconnects  permitted  the 
power  supply  leads  to  be  quickly  changed  to  either  feed-through.  Resistive  heatin< 
provided  the  necessary  heat  to  cause  evaporation  (sublimation)  of  the  coating  metal.  The 
coating  metal  was  held  in  tungsten  foil  boats  which  were  clamped  to  the  power  feed- 
throughs.  The  upper  ring  of  three  service  ports  housed  a  rotary  shutter  shield  to  cover  the 
substrate  when  applying  power  to  the  feed-throughs.  This  enabled  any  oxides  or 
hydrocarbons  to  be  burned  off  the  coating  metal  prior  to  actually  depositing  the  metal  on 
the  substrate  surface.  Once  the  burnoff  was  completed,  the  shutter  could  be  rotated  out  of 
the  way  to  expose  the  substrate  surface.  The  second  upper  service  port  housed  the  film 
thickness  monitor  coaxial  cable  feed-through.  The  remaining  upper  service  port  housed  a 
combination  electrical  feed-through  which  passed  a  K  type  thermocouple  lead  for 
monitoring  the  substrate  surface  temperature  during  metallization,  and  two  DC  power 
leads  used  to  supply  power  to  a  small  mica  heater  unit.  The  mica  heater  unit  capable  of 
heating  the  substrate  up  to  500°C  during  metallizations.  The  substrate  itself  was 
physically  clamped  to  the  heater  backing  plate  and  held  at  the  desired  position  by  a  small, 
adjustable  stand  located  inside  the  chamber.  All  seals  were  made  with  rubber  o-rings 
coated  with  high  vacuum  grease. 

The  high  vacuum  needed  to  perform  metallizations  was  achieved  by  use  of  a 
diffusion  pump  which  was  capable  of  reaching  vacuum  levels  as  low  as  5  x  10"7  torr  under 
ideal  conditions.  Continuous  pressure  readings  of  the  vacuum  during  metallizations  was 
provided  by  use  of  a  Penning  gage  with  a  digital  readout 
2.      Film  Thickness  Monitor 

The  film  thickness  monitor  provided  instantaneous  readings  of  the  film  thickness 
during  metallizations.  The  film  thickness  monitor,  supplied  by  Cressington  Inc.,  used  a 
pizo-electric  quartz  crystal  located  inside  the  chamber.  The  signal  developed  was  sent 
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through  the  coaxial  cable  to  a  processing  unit  with  a  digital  readout.  Since  the  oscillating 
characteristics  of  the  crystal  is  a  function  of  the  coating  metal  density,  the  density  of  the 
coating  metal  was  required  as  an  input  to  the  film  thickness  monitor  processor  prior  to 
deposition.  Additionally,  the  location  of  the  crystal  (distance  and  angular  separation)  with 
respect  to  the  substrate  and  the  coating  source,  greatly  affected  the  thickness  readings 
obtained.  This  geometry  effect  was  accounted  for  by  use  of  a  scaling  factor,  called  a 
tooling  factor,  which  was  also  a  required  input  to  the  thickness  monitor  processor.  Figure 
1 1  shows  the  key  dimensions  and  angles  which  required  accurate  measurement  in  order  to 
calculate  the  correct  tooling  factor.  The  film  thickness  monitor  processor  used  was  a  two 
channel  device  that  allowed  two  different  coating  metals  to  be  monitored  during 
deposition.  The  two  channel  feature  was  necessary  in  order  to  monitor  the  inter-layer 
thickness  as  well  as  the  Au  film  thickness. 

As  shown  in  Figure  1 1,  the  distances,  ls  and  ld,  greatly  influenced  the  magnitude 
of  the  tooling  factor.  Selection  of  the  correct  sample-to-source  distance,  ls ,  and  the 
crystal-to-source  distance,  ld,  was  critically  important  in  other  regards  as  well: 

(1)  Too  small  a  value  of  ls  resulted  in  splattering  of  the  coating  metal  on  the 
substrate  surface, 

(2)  Too  small  a  value  of  ld  resulted  in  the  rapid  coating  of  the  crystal  face  to  the 

point  that  the  crystal  was  rendered  useless  before  achieving  the  desired  film  thickness, 

(3)  Too  large  a  value  of  ls  resulted  in  expenditure  of  a  large  amount  of  coating 
metal,  greatly  lengthening  the  time  required  to  coat  the  substrate  which  resulted  in 
excessive  heating  of  the  chamber  and  erroneous  thickness  readings.  Unfortunately,  the 
optimum  values  of  ls  and  ld  had  to  be  obtained  through  trial  and  error.  In  addition,  the 
optimum  distances  differed  significantly  depending  on  the  coating  metal  used.  For  the 
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Figure  11.  Calculation  of  the  Tooling  Factor  From  Metallization  Geometry  [Ref.  66] 
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metals  used  in  this  study,  an  ls  of  10.7  cm  and  an  ld  of  28.0  cm  were  found  to  be 
satisfactory  for  metallization. 
C.     PROCEDURE 

1.  Gold  Standard  Preparation 

Preparation  of  an  Au  standard  was  necessary  to  obtain  the  instrumental  profile 
used  in  the  Warren- Averbach  analysis  as  discussed  previously.  The  Au  standard  was 
prepared  by  annealing  a  sufficient  amount  of  the  Au  powder  in  flowing  Argon  at  300°C 
for  three  hours  in  an  alumina  crucible.  After  air  cooling,  it  was  found  that  the  Au  particles 
had  sintered  together.  To  break  up  the  particles,  the  Au  powder  was  shaken  in  an  ultra- 
sonic vibrator  and  then  sieved  using  a  325  mesh  sieve  onto  a  low  noise,  quartz  XRD 
sample  holder.  Adhesion  of  the  powder  to  the  sample  holder  was  achieved  by  using  a  thin 
film  of  high  vacuum  grease.  A  glass  slide  was  then  lightly  pressed  over  the  powder  to 
assure  that  the  powder  was  as  flat  as  possible. 

2.  Substrate  Surface  Preparation 

Table  2  identifies  the  samples  by  number  and  provides  a  metallization  summary. 

TABLE  2.  SUMMARY  OF  SURFACE  TREATMENTS 


SAMPLE 
NAME 

BASE 

PRESS. 

torr 

SURF. 
ROUGH, 
microns 

INTER- 
LAYER 
metal/nm 

OXIDE 
LAYER 

microns 

TEMP. 

DURING 

DEP. 

degC 

BAKED 
(Prior  to 
Coating) 

FINAL 

THICK. 

nm 

AUALN1 

5xl0~6 

20-25 

NONE 

NONE 

25-60 

NO 

327 

AUALN2 

5xl0-6 

15 

NONE 

NONE 

25-60 

NO 

328 

AUALN3 

5xlCf6 

NONE 

NONE 

25-60 

NO 

328 

AUALN4 

4xl0"6 

Cr/6.4 

NONE 

25-60 

YES 

310 

AUALN5 

6xl0"6 

Ti/7.7 

NONE 

25-60 

YES 

320 

AUALN6 

4.5  xlO-6 

NONE 

NONE 

25-60 

YES 

327 

AUALN7 

4.5  xlO"6 

N(  )NE 

3.0 

25-60 

YES 

327 

AUALN8 

6xl(T6 

Cr/6.3 

noni-: 

95-140 

YES 

329 

AAU1 

5  xlO-9 

20-25 

NONE 

NONE 

25-60 

NO 

450 
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Prior  to  preparing  the  substrate  surfaces,  the  A1N  substrates  were  cut  to  a 
dimension  that  would  fit  directly  in  the  XRD  sample  holder  using  a  diamond  saw.  The 
surface  roughness  condition  of  the  AUALN1  and  AAU1  samples  were  used  "as  received" 
and  estimated  to  be  between  20  and  25  microns.  No  other  treatments  were  applied  to 
these  two  samples.  The  AUALN2  sample  was  prepared  using  first  400  then  600  grit 
abrasive  paper  with  de  ionized  water.  The  remaining  samples,  having  a  surface  roughness 
of  1  micron,  were  prepared  by  using  400  and  600  grit  abrasive  paper,  followed  by 
polishing  with  6  micron,  then  1  micron  diamond  paste.  Prior  to  metallization,  all  samples 
were  thoroughly  washed  with  acetone  followed  by  ultrasonic  sink  cleaning  in  ethanol. 

The  AUALN7  sample  was  oxidized  using  the  procedure  of  [Ref.  67].  The  A1N 
substrate  was  placed  in  a  fused  silica  holder  and  heated  to  1000°C  using  a  dilatometer 
oven.  Clean,  dry  oxygen  was  supplied  at  a  flowrate  of  20  ml/hr  (resulting  in  a  linear 
flowrate  of  about  40  mm/min)  for  a  period  of  15  hours.  A  thermocouple  was  used  to 
monitor  the  substrate  temperature  during  oxidation.  Based  on  these  conditions,  an  oxide 
layer  of  about  3  microns  was  produced  on  the  substrate  surface. 
3.      Metallization  Procedure 

After  preparing  and  cleaning  the  substrate  surfaces,  the  substrates  were  clamped 
to  a  combination  sample  holder/heater  assembly  in  the  vacuum  chamber.  Vacuum  was 
then  drawn  in  the  chamber  with  the  diffusion  pump.  An  average  vacuum  of  5  x  10"6  torr 
was  achieved  during  all  metallizations  with  the  exception  of  the  AAU 1  sample  which  was 
prepared  at  5  x  10"9  torr  in  a  separate  UHV  chamber.  A  thermocouple  attached  to  the 
sample  holder  was  placed  in  contact  with  the  edge  of  the  substrate  which  allowed  for  the 
continuous  monitoring  of  the  substrate  temperature  during  metallization.  The  rotary 
shutter  was  then  inserted,  isolating  the  substrate  from  the  metal  source.  Power  was 
supplied  to  each  power  feed-through  to  burn  off  any  oxides  or  hydrocarbons  from  the 
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coating  metal.  The  strip  heater  was  energized,  using  a  variable  DC  power  supply,  until  the 
substrate  was  heated  to  300°C  for  a  period  of  30  minutes.  This  baking  procedure  was 
used  on  all  samples  except  samples  AUALN1  through  AUALN3,  since  the  heater  unit  was 
not  available  when  these  samples  were  prepared.  The  baking  was  performed  as  a 
precaution  to  remove  any  trace  amounts  of  hydrocarbons  on  the  substrate  surface  prior  to 
coating.  With  the  exception  of  the  AUALN8  sample,  the  substrate  was  allowed  to  return 
to  a  temperature  of  25° C.  The  AUALN8  sample  was  allowed  to  cool  to  95" C  to 
investigate  the  substrate  surface  temperature  effect  on  film  quality.  The  rotary  shutter  was 
then  removed  and  the  metallization  of  the  substrate  was  performed. 

Inter-layers  of  Cr  and  Ti  (samples  AUALN4,  AUALN5  and  AUALN8)  were 
made  during  the  metallization  process  by  using  one  of  the  two  power  feed-throughs  to 
deposit  the  inter-layer  metal.  The  remaining  power  feed-through  was  then  used  to  deposit 
the  Au  film.  Since  two  power  feed-throughs  were  available,  chamber  vacuum  did  not 
have  to  be  broken  between  the  deposition  of  the  inter-layer  metal  and  the  Au  film.  This 
prevented  the  oxidation  of  the  inter-layer  surface  prior  to  applying  the  Au  film. 

During  the  inter-layer  and  Au  metallizations,  an  average  deposition  rate  was  280 
angstroms  per  minute  was  used.  A  deposition  temperature  for  samples  AAU1  and 
AUALN1  through  AUALN7  was  maintained  between  25-60°C.  Deposition  temperature 
for  the  AUALN8  sample  was  maintained  between  95  -  140°C. 

Once  metallization  was  complete,  the  rotary  shutter  was  inserted  prior  to 
breaking  chamber  vacuum,  to  prevent  blowing  debris  onto  the  freshly  coated  substrate 
surface. 
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D.     XRD  PROCEDURES 

1.  Parameters  For  XRD  Scanning 

Following  metallization,  the  samples  were  loaded  directly  into  the 
diffractometer.  A  Philips  XRG  3 100  X-ray  generator  with  a  copper  target,  a  bent  graphite 
monochromator  and  a  Digital  Systems  data  control  and  processor  were  used  to  perform 
the  scans.  Goniometer  scans  were  performed  using  a  step  size  of  0.01  degrees  (2-theta) 
for  4  seconds.  Tube  voltage  was  maintained  at  30  kv  with  a  tube  current  of  35  milliamps. 
The  same  settings  were  used  to  scan  the  Au  standard  and  all  samples. 

2.  Post  XRD  Scan  Processing 

Following  the  XRD  scans,  the  raw  data  scan  files  generated  were  analyzed  using 
a  APD  1700  software  package  developed  by  Philips.  The  capabilities  of  this  software 
enabled  the  raw  data  to  be  processed  in  a  way  suitable  for  a  Warren- Averbach  analysis 
described  earlier.  Specific  software  capabilities  discussed  below,  were  taken  from  [Ref. 
68]. 

The  first  step  of  the  analysis,  called  the  Pattern  Treatment  module,  was  used  to 
separate  crystalline  peaks  from  amorphous  peaks.  Next,  K-alpha  2  lines  and  background 
were  removed  and  a  new  file  was  generated.  At  this  point,  the  operator  had  the  option  to 
modify  the  new  file,  by  selecting  an  intensity  level  below  which  the  system  ignored  any 
intensity  information.  This  feature  enabled  the  operator  to  filter  out  any  background  noise 
spikes.  It  was  found  that  this  parameter  played  a  major  role  in  obtaining  good  profile  fits 
of  the  raw  data. 

Following  the  Pattern  Treatment  phase  of  the  analysis,  the  raw  data  files 
developed  were  used  to  generate  a  mathematical  profiles  of  the  (1 1 1)  and  (222)  Au  peaks 
occuring  in  the  coated  A1N  samples.  This  part  of  the  analysis,  called  Profile  Fitting 
module,  employed  a  Marquardt  non-linear  least  squares  algorithm  to  fit  the  profiles. 
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Profile  fitting  was  necessary  to  separate  the  peak  overlap  between  A1N  and  Au  peaks. 
Figures  12  and  13  show  diffraction  scans  of  the  Au  standard  and  A1N  prior  to 
metallization.  Note  that  both  the  (1 1 1)  and  (222)  peak  positions  for  Au  had  interfering 
AIN  peaks  that  required  separation  before  the  Warren- Averbach  analysis  could  be 
performed.  The  Profile  Fit  module  generated  profiles  for  both  the  AIN  and  the  Au  peaks 
present.  After  the  profiles  were  generated,  the  operator  could  super  impose  the  Au  and 
AIN  profiles  onto  the  raw  data  scan  as  well  as  any  residual  peaks  that  were 
"mathematically"  left  over  from  the  profile  fitting  process.  The  goodness  of  fit  was 
checked  by  observing  the  intensity  of  the  residual  peaks  in  the  2-theta  range  of  interest.  If 
the  residuals  showed  low  intensities  over  the  2-theta  range  of  interest,  (i.e.,  if  no  peaks 
other  than  the  expected  AIN  peak  and  the  expected  Au  peak  were  apparent),  the  fit  was 
considered  good.  Figures  14  and  15  show  the  result  of  this  procedure  for  an  Au  (1 1 1) 
peak  and  a  (222)  peak  respectively.  It  can  be  seen  from  Figures  14  and  15  that  only  a  AIN 
peak  and  a  Au  peak  are  present  and  no  others.  The  profile  fitting  method  was  largely  trial 
and  error,  using  the  desired  2-theta  range,  peak  intensity,  peak  position  and  background  as 
input  parameters.  Once  satisfactory  profiles  were  developed  for  both  the  Au  (111)  and 
(222)  peaks,  they  were  stored  in  separate  files  to  be  used  in  the  Warren-Averbach  phase  of 
the  analysis. 

The  Warren-Averbach  module  used  the  Au  (1 1 1)  and  (222)  profiles  of  a  sample 
as  well  as  the  (1 1 1)  and  (222)  peaks  of  the  Au  standard  and  performed  the  Warren- 
Averbach  deconvolution  using  either  the  cosine  method  or  the  modulus  method.  Both 
methods  were  used  in  this  study  and  compared.  Particle  size  and  r.m.s.  strain  are  the 
output  values.  In  this  module,  corrections  for  insu"umental  errors,  the  Lorentz- 
Polarisation  factor,  the  K  alpha-2  component  and  background  corrections  were  applied. 
Once  the  corrections  were  applied,  the  deconvolution  process  for  the  low  order  ((111) 
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peaks)  was  executed.  The  program  performed  a  Fast  Fourier  Transform  (FFT)  on  both 
peaks  followed  by  a  Stokes  deconvolution  algorithm.  The  Fourier  coefficients  for  'he  true 
profile  were  then  computed  and  stored.  The  entire  process  was  repeated  for  the  high 
order  (222)  peaks.  Once  the  Fourier  coefficients  for  both  order  reflections  were  obtained, 
the  module  performed  the  necessary  calculations  and  determined  the  particle  size  and 
r.m.s.  strain  values  for  the  sample.  A  hook  correction,  if  needed,  was  performed  using  an 
automatic  method  or  a  manual  method. 
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IV.  RESULTS  AND  DISCUSSION 

A.     RESULTS  SUMMARY 

1.      R.M.S.  Srain,  Particle  Size  and  Stress  Deviation 

Table  3  lists  the  results  of  the  panicle  size  and  r.m.s.  strain  obtained  by  the 
Warren-Averbach  analysis.  Results  are  shown  for  both  the  cosine  method  and  the 
modulus  method  and  show  fairly  good  agreement  for  all  samples.  Table  3  also  lists  the 
stress  deviation  which  was  calculated  based  on  [Ref.  69].  This  method  assumes  that  the 
stress  state  of  the  film  is  simply  biaxial,  since  the  film  thickness  is  very  small  with  respect 
to  film  length  and  width.  The  biaxial  stress  equation  used  for  this  calculation  is 


Kv)  =  — (^>r 


where:  (odev)  represents  the  possible  stress  deviation  from  the  mean  stress,  E  is  Young's 
Modulus,  o  is  Poisson's  ratio  and  (e, )       is  the  r.m.s.  strain.  Values  for  Young's  Modulus 

and  Poisson's  Ratio  were  taken  from  [Ref.  70]  and  represent  thin  film  values  vice  bulk 
values.  The  values  of  r.m.s.  strain  found  in  this  study  are  of  the  same  order  of  magnitude 
as  r.m.s.  strains  reported  in  other  studies  [Ref.  71]  for  Au  films  of  roughly  the  same 
thickness. 

2.      Discussion 

Considering  the  surface  roughness  results  (AUALN1,  AUALN2  and  AUALN3), 
it  was  observed  that  as  the  surface  roughness  of  the  substrate  was  increased,  the  r.m.s. 
strain  increased.  Since  this  series  of  films  were  made  under  the  same  metallization 
conditions  and  involve  the  same  materials,  the  (ACTE)(AT)  for  each  sample  was  the  same. 
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This  indicates  a  greater  extent  of  plastic  deformation  occurs  on  rougher  substrates  and 
may  indicate  that  adhesion  improves  as  the  substrate  surface  roughness  is  increased.  A 
likely  explanation  for  this  behavior  is  that  better  mechanical  interlocking  is  possible 
between  the  film  and  the  substrate  as  surface  roughness  increases.  This  interlocking 
reduces  the  degree  of  strain  relaxation  possible  by  interfacial  sliding,  thereby,  enhancing 
the  extent  of  plastic  strain. 

The  two  samples  metallized  under  different  vacuum  conditions  (AUALN1  and 
AAU1)  are  considered  next.  Both  samples  had  the  same  surface  roughness  and  were 
metallized  under  the  same  conditions  with  the  exception  of  the  metallization  pressure.  The 
AAU1  sample  (metallized  under  UHV  conditions)  had  an  r.m.s.  strain  significantly  higher 
than  the  AUALN1  sample  (metallized  under  HV  conditions),  indicating  that  more  plastic 
deformation  occurred  in  the  film  depositied  at  a  lower  base  pressure.  Since  the  same 
materials  were  used,  it  is  suspected  that  the  adhesion  improves  as  metallization  pressure  is 
reduced.  A  possible  explanation  for  this  behavior  is  that  less  contamination  of  the 
substrate  surface  occurs  at  lower  metallization  pressures,  resulting  in  better  interfacial 
bonding.  Other  studies  of  copper  films  have  borne  this  out  [Ref.  72].  Additionally,  lower 
metallization  pressures  are  known  to  promote  better  epitaxial  film  growth,  which  may  be 
related  to  the  degree  of  interfacial  bonding. 

Analysis  of  the  interlayer  results  is  a  bit  more  complex  since  two  interfaces  need 
to  be  considered.  If  the  interlayer  is  assumed  to  act  as  a  compliant  layer  (i.e.,  provides  no 
resistance  to  the  shrinkage  of  the  applied  film  during  cooling),  then  no  stress  could  be 
transferred  from  the  interlayer  to  the  film.  If  the  interlayer  acts  as  a  separate  entity,  then 
the  differences  in  the  CTE's  between  the  interlayer  and  the  film  should  account  for  the 
degree  of  plastic  deformation,  hence,  the  degree  of  r.m.s.  strain  observed.  Finally,  the 
interlayers  may  be  considered  to  be  merely  surface  modifiers,  since  the  thickness  of  the 
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interlayer  is  very  small  compared  to  the  film  thickness.  If  interlayers  do  act  as  surface 
modifiers,  then  the  (ACTE)(AT)  effect  between  the  interlayer  and  the  film  is  not 
considered,  and  the  presence  of  the  interlayer  either  improves  or  degrades  interfacial 
bonding  between  the  film  and  the  substrate  in  some  fashion. 

Analysis  of  the  samples  with  inter-layers  (AUALN4,  AUALN5  and  AUALN8) 
show  significant,  measurable  r.m.s.  strain  values,  indicating  that  the  inter-layers  probably 
do  not  act  as  compliant  layers.  The  Au/AIN  sample  with  a  Cr  inter-layer  (AUALN4) 
displays  a  higher  r.m.s.  sttain  than  the  Au/AIN  sample  of  the  same  surface  roughness 
without  any  inter-layer  (AUALN6),  although  the  CTE  mismatch  between  Au  and  Cr  is 
less  than  that  between  Au  and  A1N.  Therefore,  it  is  unlikely  that  the  inter-layers  act  as 
separate  entities.  This  is  reasonable  since  the  inter-layers  used  were  very  thing, 
-60-70  A ,  and  would  not  be  expected  to  act  as  substrates.  It  is  believed  that  the  inter- 
layers act  as  surface  modifiers,  either  by  altering  the  nature  of  chemical  bonding  between 
Au  and  A1N  at  the  interface,  or  by  forming  good  chemical  bonds  with  both  the  film  and 
the  substrate  materials.  The  results  show  that  the  r.m.s.  strain  in  the  Au  film  is 
significantly  larger  for  the  sample  with  the  Cr  inter-layer  than  for  the  sample  without  any 
inter-layer,  indicating  that  Cr  promotes  interfacial  adhesion.  The  exact  reasons  for  this  is 
not  known  and  would  require  further  study.  It  is  suspected  that  Cr  has  a  high  affinity 
toward  oxygen,  and  any  oxygen  which  may  have  been  absorbed  on  the  surface  of  the  A1N 
could  result  in  good  bonding  between  Cr  and  A1N  whereupon  strong  metallic  bonds  could 
be  formed  between  Cr  and  the  Au  film.  Likewise,  it  has  been  reported  that  Ti  forms 
nitrides  at  the  interface  with  A1N,  allowing  the  Ti  to  bond  strongly  with  the  substrate  [Ref. 
73].  When  a  metal  is  then  deposited  on  Ti,  metallic  bonds  are  formed  between  the  film 
and  the  Ti  inter-layer,  giving  rise  to  a  film  well  bonded  to  the  substrate.  Accordingly,  one 
would  expect  that  the  r.m.s.  strain  measured  in  the  sample  with  the  Ti  inter-layer  would  be 
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larger  than  that  in  the  sample  without  an  inter-layer  (AUALN3).  However,  as  seen  from 
Table  3,  the  sample  with  a  Ti  inter-layer  had  a  r.m.s.  strain  much  lower  than  the  others. 
This  result  may  be  due  to  oxygen  contamination  of  the  Ti  inter-layer  during  deposition, 
thereby  preventing  strong  metallic  bonding  between  the  Au  and  A1N  and  resulting  in  low 
levels  of  plastic  strain.  Further  chemical  investigation  of  the  interface  regions  are  required 
to  verify  this  hypothesis. 

The  Cr  inter-layer  sample  that  was  metallized  at  a  higher  substrate  temperature 
(AUALN8)  shows  a  much  lower  degree  of  r.m.s.  sttain  than  for  the  Cr  inter-layer  sample 
(AUALN4)  which  was  metallized  at  a  lower  substrate  temperature.  A  possible 
explanation  for  this  observation  is  that  the  inter-layer  of  the  sample  metallized  at  the 
higher  substrate  temperature  could  achieve  a  greater  extent  of  stress  relaxation  (through 
increased  dislocation  movement)  prior  to  the  application  of  the  Au  film.  This  could  result 
in  less  plastic  stress  being  transferred  to  the  Au  film  and  less  r.m.s.  strain  upon  cooling. 

Finally,  the  substrate  treated  with  an  oxide  layer  (AUALN7),  showed  a  significant 
amount  of  r.m.s.  strain.  A  possible  explanation  is  that  Au-O  bonding  occurred  at  the 
oxide-film  interface  when  the  high  temperature  Au  metal  first  reached  the  oxide  surface 
during  metallization.  Since  it  is  known  that  oxygen  diffuses  extremely  well  into  A1N,  good 
bonding  at  both  interfaces  is  possible,  resulting  in  a  large  degree  of  plastic  stress  transfer 
to  the  Au  film  and  high  r.m.s.  strain.  This  hypothesis  would  also  require  further  chemical 
analysis. 

The  average  particle  size  obtained  from  the  Warren- Averbach  analysis  is  a  function  of 
dislocation  spacing,  and  the  degree  stacking  faults  and  twinning.  If  the  degree  of  stacking 
faults  and  twinning  are  the  same  for  all  the  samples,  then  the  particle  sizes  obtained  in  this 
study  are  indicative  of  the  subgrain  size  of  Au  film.  Obtaining  stacking  fault  and  twinning 
probabilities  from  X-ray  diffraction,  through  the  measurement  of  peak  shifts  relative  to  a 
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strain  free  standard  of  sufficient  particle  size,  is  possible  [Ref.  74],  and  further  studies 
using  these  methods  would  provide  more  information  on  the  microstructure  of  the  Au 
films  produced. 

3.      Uncertainty 

The  APD  1700  software  used  to  perform  the  Warren- A verbach  analysis 
performs  uncertainty  estimates  as  well.  Table  4  lists  the  uncertainties  of  both  r.m.s.  strain 
and  particle  size  for  each  method  used.  It  should  be  noted  that  the  uncertainties  listed 
apply  to  the  computational  methods  used  in  the  Warren-Averbach  analysis  phase  of  the 
APD  1700  program  only,  and  do  not  represent  the  combined  uncertainties  of  the 
metallization  process  or  from  the  XRD  measurements. 
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TABLE  4.  UNCERTAINTY  SUMMARY 


SAMPLE 

COSINE 

METHOD 

MODULUS 

METHOD 

R.M.S.  STRAIN 

% 

PARTICLE  SIZE 

% 

R.M.S.  STRAIN 

% 

PARTICLE  SIZE 

% 

AUALN1 

3.5 

4.2 

3.0 

3.9 

AUALN2 

19.2 

14.7 

2.96 

3.01 

AUALN3 

30.8 

15.6 

23.1 

14.5 

AUALN4 

27.6 

21.1 

19.4 

14.5 

AUALN5 

4.0 

3.0 

4.4 

4.0 

AUALN6 

24.5 

18.0 

24.7 

17.6 

AUALN7 

34.1 

15.5 

33.6 

15.0 

AUALN8 

40.0 

24.6 

13.7 

S.16 

AAU1 

17.9 

18.2 

18.7 

17.8 
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V.  CONCLUSIONS 

The  r.m.s.  strain  value  increased  as  the  surface  roughness  of  the  substrates  was 
increased  indicating  a  higher  degree  of  plastic  deformation  for  roughened  substrates.  It  is 
suspected  that  a  higher  surface  roughness  improves  mechanical  interlocking  between  the 
film  and  the  substrate  which  improves  adhesion. 

A  lower  metallization  pressure  resulted  in  a  much  higher  value  for  r.m.s.  strain,  which 
indicates  better  adhesion  is  possible  by  improving  vacuum  during  deposition.  Less 
contamination  of  the  substrate  surface  is  a  possible  reason  for  this  observation. 
Additionally,  deposition  pressure  may  effect  the  epitaxial  growth  of  the  film  crystals  which 
may  be  related  to  film  adhesion. 

The  presence  of  an  oxide  layer  was  seen  to  raise  the  r.m.s.  strain  values  significantly, 
indicating  the  possibility  of  metal-0  bonding  between  the  Au  and  the  Alumina  layer. 
Although  improved  adhesion  is  indicated  with  the  presence  of  an  oxide  layer,  the  thermal 
conductivity  of  the  A1N  /  A1203  /  Au  system  is  most  likely  lower  than  the  A1N  /  Au  system 
due  to  oxygen  diffusion  into  the  A1N. 

Measurable  r.m.s.  strains  for  metallizations  using  inter-layers  indicated  that  the  inter- 
layer  does  not  act  as  a  compliant  layer.  Additionally,  the  inter-layer  is  not  likely  to  be 
acting  as  a  separate  entity  since  the  ACTE  between  the  inter-layer  and  the  Au  film  does  not 
account  for  the  relative  amounts  of  r.m.s.  strain  observed  for  the  Cr  and  Ti  inter-layer 
samples.  It  is  believed  that  inter-layers  act  as  surface  modifiers  that  change  the  bonding 
characteristics  of  the  interface  in  some  way.  A  high  r.m.s.  strain  value  was  observed  using 
a  Cr  inter-layer,  indicating  it  improves  adhesion  compared  to  a  system  with  no  inter-layer. 
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A  Ti  inter-layer  was  seen  to  produce  a  small  r.m.s.  strain  value  indicating  less  adhesion, 
possibly  due  to  contamination  of  the  Ti  by  oxygen  diffusion  during  deposition. 
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VI.  RECOMMENDATIONS 

The  following  recommendations  are  given  for  further  study  of  metallized  A1N 
substrates: 

•  More  Au  thin  films  deposited  on  A1N  substrates  using  the  same  substrate  preparation 
parameters  and  metallizing  conditions  should  be  performed  to  improve  the  statistical 
reliability  of  these  results. 

•Careful  analysis  of  the  stacking  fault  and  twinning  probabilities  from  the  X-ray  data, 
and  the  extraction  of  real  crystallite  sizes  should  be  used  to  supplement  the  r.m.s. 
strain  data. 

•  TEM  studies  of  the  prepared  films  would  be  extremely  useful  in  determining  the 
amount  of  stacking  and  twin  faults  which  could  be  used  to  verify  the  effective  particle 
size  results  obtained  in  this  study. 

•  The  films  prepared  should  be  checked  for  adhesive  strength  by  scratch  or  peel  testing 
to  allow  the  r.m.s.  strain  values  obtained  in  this  study  to  be  directly  correlated  to 
adhesive  strength.  If  a  direct  correlation  indeed  exists,  then  the  r.m.s.  strain  values 
could  be  used  directly  to  predict  thin  film  performance  without  the  need  for 
destructive  testing. 

•  Post-metallization  heat  treatments  may  hold  promise  in  reducing  residual  stress 
levels,  through  dislocation  movement,  as  well  as  allowing  metal  diffusion  between  the 
substrate,  interlayer  and  film.  Metal-to-metal  bonding  may  also  be  improved  at 
elevated  metallization  temperatures  by  increased  diffusion.  A  systematic  study  of 
substrate  metallization  temperature  and  post-metallization  heat  treatments  is  therefore 
warranted. 
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